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ABSTRACT 

Let  z(t)  €  r”  be  a  generalized  Poisson  process  with  parameter  X. 

In  the  present  paper,  the  conditions  of  existence  and  limiting  behavior 
as  X-*-*  or  as  X-»0  of  the  stationary  distribution  of  the  solution 
of  Langevin  equation  dx(t)  ■  Ax{t)  +  dz(t)  are  investigated.  Using  these 
results,  the  distribution  of  virtual  waiting  time  in  a  queueing  system  with 
variable  service  speed  is  studied. 
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1.  Introduction 

Let  z(t)  =  (z, (t),  z^(t)...,  z  (t))^  e  be  a  generalized 

X  ^  II 

Poisson  process  with  parameter  X  and  jun^s  x^,  X2,...,  Xjjj,....  Let  also 
A:  -►  be  a  linear  onerator  determined  by  the  matrix  A  =  ||  a- •  ||  n 

„  •  i»i=l 

(We  siqjpose  that  the  basis  in  k  is  fixed.) 

The  present  paper  deals  with  conditions  of  existence  and  limiting  pro¬ 
perties  as  X  “  or  as  X  -►  0  ,  of  stationary  distribution  of  the  process 
x(t)  =  (Xj^(t),...,  x^(t))^  e  ,  which  satisfies  the  formal  equation 

dx(t)  =  Ax(t)dt  +  dz(t)  (1.1) 

The  equation  (1.1)  is  a  continuous  analog  of  the  autoregression  eauation 

Vl  ■  *  Vl 

where  the  distribution  of  independent  identically  distributed  (i.i.d.) 
random  vectors  z ^ ,  ^2’"’  weight  1  -  e  .  Thus , 

it  can  describe  some  processes  connected  with  radioactive  decay,  queuing 
systems  with  changeable  rate  of  service,  etc. 

2.  Stationary  Distribution  and  the  Conditions  of  its  Existence 

As  usual,  we  will  assume  that  the  differential  equation  (1.1)  with 
initial  condition  x(0)  =  Xq  is  equivalent  to  the  integral  equation. 

x(t)  “  Xg  +  j  Ax(u)du  +  z(t)  (2.1) 

0 

vdiich  holds  with  probability  one  for  all  values  of  t. 

In  what  follows,  we  assume  that  the  process  z(t)  has  rightcontinuous 
sanple  paths  with  probability  one. 


Lenina  2.1.  The  eqioation  (2.1)  has  a  unique  solution  in  the  class  of 
measurable  processes.  This  solution  is  a  rightcontinuous  strongly  Markovian 
process  and  can  be  written  in  the  form 

x(t)  =  exp  {At}xQ  +  f  exp  { A(t-u) }dz(u)  (2.2) 

0 

where  the  integral  on  the  right-hand  side  of  (2.2)  is  a  Stilt jes  integral 
and  exists  with  probability  one. 

The  proof  of  this  lemna  is  routine. 

Lemna  2.2.  The  one-dimensional  distributions  of  the  process  x(t) 
are  infinitely  divisible  and  have  the  characteristic  functions  (c.f.) 

4'(s;t)  =  E  exp  {i(s,x(t))}  = 

exp  {i(s,exp  {At}xQ)  -  xj  (exp  {A^u}s))du}  (2.3) 

0 

where  (s)  *  E  (exp  i(s,Xj^)}  . 

The  proof  of  this  lemma  follows  from  the  representation  (2.2). 

Theorem  2.1.  The  process  x(t)  possesses  limiting  distribution  as 
t  -►  00  vdiich  does  not  depend  on  the  initial  state  x^  if  and  only  if 

1)  the  eigenvalues  of  A  lie  in  the  left  halfplane; 

2)  E  log  (1+lXj^l)  <  00  . 

Proof.  If  ’l'(s;t)  =  'l'(s;t,XQ)  -►  4'(s)  as  t  -►  oo  and  'l'(s)  is  continuous, 
then  (s)  does  not  vanish  since  it  is  a  c.f.  of  an  infinitely  divisible 
distribution  in  .  Thus 

ejq)  {i(s,exp  {At}xQ)  »  ’i'(s;t,2xQ)T’^(s;t,XQ)  -*•  1 

oo 

for  all  initial  values  Xg  and  s  €  1^  . 

This  is  possible  if  and  only  if  the  condition  1)  holds.  In  this  case. 


we  have  the  equality 
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I  T 

4’(s)  =  exp  {-X  j  fl'i^fexp  {A^u}s))du}  (2.4) 

0 

It  is  easy  to  check  that  the  infinitely  divisible  c.£.  'i'(s;t)  and 

4'(s)  deteimined  by  (2.3)  and  (2.4)  have  the  Lfevy  representations 

log  'l'(s;t)  =  i(s,Y^)  -  Q^(s)  +  j  [exp  {i(s,x)} 

|x|>0 

-  1  -  i(s,x)(l+(x,x))’^]N^(dx)  (2.5) 

log  4'(s)  i(s,Y)  -  Q(s)  +  j  [exp  {i(s,x)} 

|x|>0 

-  1  -  i(s,x)(l+(x,x))'^]N(dx)  (2.6) 

id\ere 

=  A  /  /  exp  {Au}x(l+(exp  {Au}x,exp  {Au}x) ) ’^P{XTedx}du  +  exp  {At}Xn  , 

^  0  ^  1  ^ 

Y  =  yU)  =  a  I  I  exp  {Au}x(l+(exp  {Au}x,exp  {Au}x))'^P{x,€dx}du  , 

0  m 

Qt(s)  =  Q(s)  =  0  , 

N^(B)  =  X  /  P{exp  {Au}x^eB}du  , 

0 

N(B)  =  N(B;X)  =  x/  P{exp  {Au}x,eB}du  . 

0 

It  follows  from  the  theorem  proved  in  [4,  p.  188]  that  ^(s.t)  -*•  4'(s) 
if  and  only  if 

N^(B)  N(B)  as  t  -►  «  for  continuity  sets  B  of  N 

lying  in  ■  (x:  |x|£e}  ; 

b)  Yt  Y  as  t  -►  00  ; 
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c)  lim  lim  f  (x,x)N^(dx)  =  0  . 

0<1X|<£ 

Using  the  sciieme  of  the  proof  of  the  theorem  1  in  [5]  we  can  show  that 
the  condition  a)  is  equivalent  to  the  condition  2)  of  the  theorem. 

Since  N^(B)  £  N(B)  ,  we  have 

0<  J  |x|N^(dx)  <  /  |x|N(dx)  (2.7) 

.  0<lxl<e  0<|xl^e 

Thus,  using  the  estimations 


[  xN(dx)  =  rd^CSp  =  f  N(S^/spdr  < 

0<lxl<£  ^  0 

J  N(I^/Spdr  =  l/  /  P{lexp  {Au}x.  | >r}dudr  £ 

0  0  0 

£  00 

xj  I  P{|xJ+l  >  cr  exp  {au}}dudr  = 

0  0 

c 

/  P{a*^log  [(|xj^|+l)(rc)‘^]  >u}dudr  = 

0  0 

X  f  Ea'^log  [(lxj^l+l)(rc)"^]dr  = 

0 

Xa  ^(eE  log  (|xj^|+l)  -  c  ^  e  log  e)  (2.8) 


vdiich  are  valid  for  small  values  of  e  and  some  a  >  0  c  >  0  ,  we  can 
easily  deduce  c) .  The  condition  b)  can  be  checked  in  the  same  manner. 


3.  Limiting  Behavior  of  the  Stationary  Distribution  as  X  -*•  «» 

Although  the  formula  (2.4)  gives  the  explicit  form  of  c.f.  of  sta¬ 
tionary  distribution  of  x(t)  ,  it  is  of  interest  to  investigate  its  limiting 
behavior  as  X  -►  »  .  In  this  part  of  the  paper,  we  will  study  the  limiting 
distributions  of  random  vectors  b*^(X) (x(t)-a(X))  under  the  assunption 
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that  x(t)  has  a  stationary  distribution.  (Here  a(X)  S  and 
b(X)  >  0  are  nonrandom  functions.)  The  class  of  nondegenerate  limiting 
distributions  for  such  vectors  coincides  with  the  class  of  stable  distribu¬ 
tions  in  ,  which  were  investigated  by  P.  Levy  [3] ,  F.  Feldhaim  [1] 
and  E.  L.  Rvaceva  [4] .  It  has  been  shown  there  that  stable  distributions  in 
r"  are  infinitely  divisible  and  their  c.f.  have  the  form 

{exp  {-|s|“[Cj^(s/|sl)+iC2(s/|sl)]  +  i(6,s)}  ,0<a<2,a/l  , 

exp  {-|s|  [Cj^(s/|s|)+iC2Cs)]  +  i(e,s)}  ,  a  =  1  (3.1) 

where  Cj^(s/|s|)  =  c  J  |cos(s^)  l°‘dH(w) , 

C2(s/is|)  =  -c  tan(2*^7ro)  Jcos(s^w)  |cos(s^)  ^dH(w)  , 

C2(s)  =  2Tr'^Jcos(^)  log  (js(  lcos(s^)  ()dH(w)  , 

Be  R^  is  a  constant  vector,  w  denotes  a  point  on  the  unit  sphere  (and 
the  vector  joining  the  origin  to  it),  H  is  a  finite  measure  on  the  unit 
sphere,  and  the  domain  of  the  integration  is  the  entire  surface  of  the 
unit  sphere.  The  number  a  is  called  the  characteristic  exponent  of  the 
distribution.  If  a  =  2  ,  we  have  the  multidimensional  normal  distribution. 
Rvaceva  [4,  p.  192]  showed  that  for  the  nondegenerate  stable  laws  in 
the  L&vy  representations  of  their  c.f. 

(s)  =*  exp  {i(B,s)-2‘^Q(s)+  f  (exp  {i(s,x)}  -  1  -  i(s,x)  (l+(x,x)) ■^)dNQ(x) 

|x|>0 

have  such  characteristics: 


a)  for  a  *  2  ,  Nq(B)  is  constant,  Q(s)  =  2(s,s)Cj^(s/|s|)  (3.2) 
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b)  for  0  <  a  <  2  ,  Nq(B)  =  R'“H(W)  ,  Q(s)  =  0  (3.3) 

for  every  set  B  of  the  form  {x:  |x|  >  R  ,  w  e  W}  ,  W  being  a  subset  of 
the  surface  of  the  unit  sphere. 

Theorem  3.1.  If  for  some  suitably  chosen  nonrandom  functions, 
a(X)  e  and  b(A)  >  0  the  distribution  of  the  vector  b*^(X) (x(t) -a(X)) 
weakly  converges  as  X  «  to  a  nondegenerate  distribution  n,  then,  n  is 
a  stable  distribution  in  with  characteristic  exponent  a  ,  0  <  a  £  2  , 
and  b(X)  is  a  regularly  varying  function  with  exponent  a’^  . 

Proof.  It  follows  from  the  formula  (2.6)  that  the  c.f.  of  the  vector 
x(t)  has  the  form  fCs)  =  exp  {AK(s)}  ,  where  K(s)  does  not  depend  on 
X.  Thus,  we  can  consider  x(t)  =  x(t;X)  as  the  value  of  a  homogeneous 
process  with  independent  increments  at  the  moment  X.  This  inplies  the 
statement  of  the  theorem.  Later,  we  will  need  the  following  result. 

Lenina  3.1.  If  the  i.i.d.  vectors  Xj^  belong  to  the  domain  of  attrac¬ 
tion  of  a  stable  law  in  with  characteristic  exponent  a  then  |xj^| 
belong  to  the  domain  of  attraction  of  a  stable  law  in  r”  with  the  same 
exponent  a  ,  0  <  o  _<  2  . 

The  proof  of  this  lemma  follows  immediately  from  the  theorems  4.1 
and  4.2  of  the  work  [4]. 

Remark  3.1.  We  will  also  use  the  fact  that  the  norming  functions 
bj^(n)  and  b2(n)  for  which  the  sequences  b^^(n)  (Xj^+. .  .+Xj^-a^)  and 
b2^(n)(|xj^|  +  ...  +  |Xjjl-a^)  weakly  converge  can  be  chosen  equal, 
bj^(n)  =  b2(n)  =  b(n)  .  This  follows  from  theorem  2.3  of  the  work  [4 j . 

We  will  consider  the  cases  a  <  2  and  a  =  2  separately. 

Theorem  3.2.  The  distribution  of  the  vector  b  ^(X) (x(t;X) -a(X)) 
weakly  converges  to  the  stable  law  in  with  characteristic  exponent 
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a  ,  0  <  a  <  2  and  spectral  Levy  measure  Nq  if  and  only  if  the  distribution 
of  the  vector  belongs  to  the  domain  of  attraction  of  the  stable  law 
with  the  same  exponent  a  and  spectral  Levy  measure  M.  The  measures  Nq  and 
M  determine  each  other  uniquely  by  the  equality 

NqCB)  =  f  M(exp  {-Au}B)du  . 

Proof.  Necessity.  From  (2.6)  we  have 

log  E  exp  {i(s,b'^(A)(x(t;A)-a{A)))}  = 

iCs.a^CA))  +  J  (exp  {i(s,x)}  -  1  -  i(s,x)(l+(x,x)) ■^)dN(b(A)x)  (3.4) 

|x|>0 

where  aj(A)  =  b*^(A) (y(A) -a(A)-  J  x(l+(x,x)) '^dN(x)) 

lxl>0 

+  f  x(l+(x,x))‘^dN(b(A)x)  . 

(x(>0 

(Note  that  the  integrals  converge  because  of  inequality  (2.8).)  It  follows 
from  theorem  1.2  [4,  p.  188]  that  the  required  convergence  is  possible  if 
and  only  if 


a)  N(b(A)B)  =  N(A;b(A)B)- Nq(B) 

A-coo 


(3.5) 


where  Nq  is  determined  by  (3.3)  and  B  C  I^/S^  is  its  continuity  set; 


b)  lim  lim  f  (x,x)dN(xb(A))  =  0; 
|xl<e 


(3.6) 


c)  lim  y(A)  =  Y 
X->00 


0 


The  condition  (3,5)  inplies  the  weak  conpactness  of  the  family  of  measures 


AP{exp  {Ar}Xj^  c  b(A)B} 


(3.7) 
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in  rP/S^  ,  e  >  0  .  To  prove  this,  we  consider  the  Borel  sets  of  the  form 
B  =  {-Av}S  ,  where  S  is  a  hypersurface  in  r”  and  the  sets  exp  {-Av}S 

do  not  intersect  for  different  values  of  v.  For  such  sets  we  have 

N(b(>,)B)  =  A  j  P{exp  {Au}xj^e  bCX)B}du  = 

0 

A  J  P{x^e  b(X)  U  exp  {-Av}S}du  =  X  f  f  P{x,  e  b(X)d  exp  {-Av}S}du  = 

0  ^  u  - 

\  f  f  duP{x  e  b(X)d  exp  {-Av}S}  =  X  /  vP{x-,  e  b(X)d  exp  {-Av}S}  -*•  N^CB)  . 

0  0^  0 

00 

Similarly,  N(b(X)  exp  {-Ar}B)  =  X  j”  P{exp  {Au}x^eb(X)  exp  {-Ar}B}du  = 

0 

r"  f" 

X  j  P{x, e  b(X)  U  exp  {-Av}S}du  =  X  /  /  P{x, e  b(X)d  exp  {-Av}S}du  = 

0  vj^u+r  0  u+r  ^ 

^v-r  « 

X  j  j  duP{xj^€  b(X)d  exp  {-Av}S}  =  X  J  (v-r)P{Xj^e  b(X)d  exp  {-Av}S}  = 
r  0  r 

00 

X  J  vP{Xj^eb(X)d  exp  {-Av}S}  -  XrP{xj^€b(X)  exp  {-Ar}B}  N^Cexp  {-Ar}B)  . 
r 


Since 


CO  00 

j  vP{x^€  b(X)d  exp  {-Av}S}  _<  j  vP{x^e  b(X)d  exp  {'Av}S}  , 


the  last  two  relations  imply  the  weak  conpactness  of  the  family  (3.7)  in 
jf/S^  ,  e  >  0  . 

In  this  case  the  family  of  measures 


XP{Xj^  €  b(X)B} 

is  also  weakly  conpact  in  R^/S^.  ,  e  >  0  . 


(3.8) 
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Since  the  estimate 

00  00 

X  j  P{exp  {Au}Xj^e  b(X)B}du  =  xj  P{exp  {Au}x^eb(X)  exp  {-At}B}du  •+ 
t  0 

NpCexp  {-At}B)  <  6 

holds  for  large  values  of  t  if  B  C  I^/S^  ,  we  can  choose  from  (3.8)  a 
weakly  convergent  subsequence 

Xj^P{x^e  b(Xj^)B}  M(B)  (3.9) 

and  interchange  the  signs  of  integral  and  limit  in  the  relation 

Nq(B)  =  lim  N(b(X.  )B)  =  lim  Xj^J  P{x,  e  t)(X, )  exp  {-Au}B}du 

k»-ao  k-^«>  0 

to  obtain  the  equality 


Ng(B)  =  /  M(exp  {-Au}B)du  . 
0 


(3.10) 


Since  the  equality  (3.10)  inplies 


M(B)  =  lim  r‘^(N„(B)-Nn(exp  {-Ar}B)) 

iM-0  ^  ^ 


the  measure  M  is  determined  uniquely  and  the  sequence  (3.8)  weakly  converges 


in  R"/S^  ,  e  >  0  , 


XP{x,  e  b(X)B}  M(B)  ,  (3.11) 

X^“ 

B  C  r"/S^  is  a  continuity  set  of  M. 

Using  the  condition  (3.6)  we  can  easily  deduce  that  for  almost  all  values 
of  u 


lim 

o+O 


lim  X 
X+oo 


f  (x,x)P{exp  {Au}Xj e  b(X)dx}  =  0 


X  <e 
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This  is  equivalent  to 


lim 

e^O 


lim 

X-HX, 


I 


X  <e 


(x,x)P{xj^e  b(A)dx} 


0 


(3.12) 


Now  we  can  use  the  obvious  continuous  analog  of  theorem  2.3  in  the  work  [4] 
which  states  that  the  conditions  (3.11)  and  (3.12)  are  sufficient  for  x^^ 
to  belong  to  the  domain  of  attraction  of  a  stable  law.  Since  norming 
function  b(A)  did  not  change,  this  stable  law  has  the  same  characteristic 
exponent  a. 

Sufficiency.  Let  the  sequence  b  ^(n) (x^+. . .+Xj^-a^)  be  weakly  con¬ 
vergent  to  a  stable  law  in  with  characteristic  exponent  a  ,  0  <  a  <  2  . 
It  follows  from  lemma  3.1,  remark  3.1  and  theorem  1  [2,  p.  313]  that  in  this 
case 


P{|xj^|>x>  =  x'“L(x)  , 


where  L(x)  is  a  slowly  varying  function.  Without  loss  of  generality  b(n) 
can  be  chosen  monotone  and  satisfying  the  relation 


nL(b(n))b  “(n)  -►1  as  n  . 

Therefore,  as  it  follows  from  the  properties  of  the  regularly  varying 
functions  (see  [2,  ch.  VIII]) 

CO  00 

X  J  P{|exp  {Au}x^|>eb(A)}du  _<  X  J  P{  jxj^  |>ceb(X)  exp  {au}}du  = 
t  t 

/oo 

P{|x, |>z}z  ^dz  <  c,  Xa'^P{ |x, |>ceb(X)  exp  (at}}  = 
ceoujexpiat}  ^  ~  ^  ^ 

C2Xa'^(ceb(X)  exp  {at}) ’°‘L(c£b(X)  exp  {at})  £  C2e  ^xp  {-aat} 
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£or  sufficiently  large  values  of  X. 
Thus,  the  condition 


nP{Xj^€  b(n)B}  ^  M(B) 


(3.13) 


which  must  be  satisfied  (see  [4,  th.  2.2])  implies 

N(b  (X)B)  =  xj  P{exp  {Au}x.  e  b(X)B}du -►  J  M(exp  {-Au}B)du  (3.14) 
0  0 

To  conplete  the  proof,  we  have  to  show  that  the  condition 


lim  lim  nb*^(n)  f  (x,x)P{Xj^ e  dx)  =  0 
|x|<eb(n) 


(3.15) 


iii5)lies 


//• 

J  (x,x)P{exp  {Au}Xj^edx}  =  0  (3.16) 

^  lxl<eb(X)  0 


In  accordance  with  theorem  2.3  of  the  work  [4],  this  conclusion  and  the 
relation  3.14)  will  be  sufficient  for  our  aim.  Using  again  the  properties 
of  regularly  varying  functions,  we  have 
00 

Xb*^(X)  J  J  (x,x)P{exp  {Au}x2  e  dx}du  = 

|x|<£b(X)  0 

.7  reb(X)  ,«> 

2Xb  ^(X)  j  J  [P{|exp  {Au}Xj^|^}  -  P{|exp  {Au}Xj^l>Eb(X) }]ydudy  ^ 

0  0 

/■eb(X) 

2Xb  '^(X)  j  j  P{|exp  {Au}Xj|^y>ydudy  ^ 

0  0 

-7 

2Xb  (X)  J  j  PljXj^l  ^cyexp  {au}}ydydu  = 

0  0 

-7  -7  /■CEb(X)exp(au} 

2c"'^Xb"'^(X)  j  J  >v}v  exp  {-2au}dvdu  _< 

0  0 
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Cj^Xb'^(X)  j  (ceb(X)  ejqj  {au})^'“L(ceb(X)  exp  {au})exp  {-2au}(iu  * 
0 

fCO  ^ 

c-Xe^  j  z'^'“L(z)dz  <  C3Xe^(ecb(X))‘“L(ecb(X))  < 
ecb(X) 


for  sufficiently  large  X  and  some  >  0  ,  a  >  0  . 

The  last  inequality  enables  us  to  prove  (3.16).  The  theorem  is  proved. 

Theorem  3.3.  The  distribution  of  the  vector  b'^(X) (x(t,X) -aCX)) 
weakly  converges  to  the  normal  law  in  with  c.f.  exp  {-2”^Qq(s)}  if 
and  only  if  the  distribution  of  the  vector  Xj^  belongs  to  the  domain  of 
attraction  of  the  normal  law  with  c.f.  exp  {-2  ^Q(s)}  .  The  quadratic 
forms  Qg(s)  and  Q(s)  determine  each  other  uniquely  by  the  equality 

-00 

QgCs)  =j  Q(exp  {AV}s)du 
0 

Proof.  In  accordance  with  theorems  1.2  and  2.3  of  the  work  [4]  and 
the  formulae  obtained  above,  we  have  to  prove  that  the  relations 


N(X;b(X)B)  -  0  ,  B  c  r"/S^  , 

X-H»  ^ 


(3.17) 


lim  lim  j  (s,x)^dN(xb(X))  =  lim  lim  j  (s,x)^dN(xb(x))  =  Qq(s)  (3.18) 
£-K)  X-«»>  „  I  I.  e^O  X-«»  I. 


0<  X  <e 


0<  X  <e 


are  equivalent  to  the  relations  XP{|Xi|  >eb(X)}  -►  0 


(3.19) 


lim  lira  X  f  (s,x)^P{x,  e  b(X)dx}  =  lim  X  f  (s,x)^P{x,  e  b(X)dx}  = 

0<|x|<e  ^  0<|x|<e 


=  Q(s) 


(3.20) 


Necessity.  The  condition  (3.19)  follows  from  (3.17)  in  the  same  manner 


as  the  condition  (3.11)  followed  from  (3.5).  The  condition  (3.18)  iuplies 


the  equicontinuity  as  \  and  as  e  -►  0  with  respect  to  s  of  the 
expressions 

X  j  (s,x)^P{exp  {Au}Xj^e  b(X)dx} 

0<|x|<e 

for  aLnost  all  ue  . 

Consequently,  we  obtain  the  equicontinuity  with  respect  to  s  of  the 
expressions 

X  f  (s,x)^P{Xje  bCX)dx}  (3.21) 

0<|x|<e 

and  can  suppose  that  a  subsequence  of  (3.21)  converges  to  the  quadratic 
form  Q(s)  .  Since 
00 

\j  j  (s,x)^P{exp  {Au}Xj^e  b(x)dx}du  = 

^  0<|x(<e 

X  r  j  (exp  {A^t}s,z)^P{exp  {Au}Xj^e  b(X)dz}du -*• 

0  0<  |exp{At}z  |<e 

T 

QQ(exp  (A  t}s)  <  c  exp  {-at}(s,s)  , 
we  obtain  from  (3.18)  that 

00 

QqCs)  */  Q(exp  {A^u}s)du  (3.22) 

0 

Since  the  equality  (3.22)  inplies 

Q(s)  =  lim  r'^(Q^(s)-Q^(exp  (A^rls) 

the  sequence  (3.21)  converges  to  Q(s)  and  the  condition  (3.20)  is  fulfilled. 

Sufficiency.  If  the  conditions  (3.19)  and  (3.20)  are  satisfied,  then  in 
accordance  with  lemna  3.1  and  remark  3.1,  Xj^  belongs  to  the  domain  of 
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attraction  of  a  normal  law  in  and  b(A)  ■  ,  where  L  is  a 

slowly  varying  function.  Then  the  inverse  function  b’^(A)  has  the  repre- 
-1  2 

sentation  b  (A)  »  A  Lj^(A)  ,  where  Lj(A)  is  a  slowly  varying  function 
too.  In  this  case,  denoting  b(A)  by  u  we  obtain 

/ao 

P{|exp  {AulXj^i  >b(A)e}  < 

0 

00  00 
A  J  >ceb(A)  e.xp  {au}}du  -  a  ^b'^tu)  f  Pilx,  |  >  v}v  ^dv  ^ 

0  tep 

00 

a'^b’Hu)6  j  (vb‘^(v))'^dv  <  Cjb‘^(u)6(u^e^Lj(ue))'^  - 

Cj^6e  ^  as  !-►“>. 

Since  6  can  be  chosen  small,  N(A;b(A) (I^/S^))  •♦0  as  A  -►  ®  and  the 
condition  (3.17)  is  satisfied.  We  divide  the  rest  of  the  proof  into  two 
parts . 

1.  E(Xj^,X2^)  <  ao  .  In  tl’is  case,  E(x(t,A)  ,x(t,A))  <  oo  and  conse¬ 

quently  x(t,A)  belongs  to  the  domain  of  attraction  of  a  normal  law.  Indeed, 
we  have  to  check  that 

J  f  (x,x)P{exp  {Au}x,  e  dx}du  <  "»  . 

0  Rfi 


We  have 
00 

/  /  (x,x)P{e3q)  {Au}Xj«  dx}du  *  2 


0  Rn 

zff 


P{ [Xj I  >  cy  exp  (au) lydydu 


rr 

0  0 


0  0 

2  J  ^  P(|xj|  >v}vc'^exp  {-2au}dvdu  «  c 


Pljexp  {Au}Xj|  >y}ydydu 


^(2a)'^E(x^,.Xj)  <  »  . 
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2.  E(Xj^,Xj^)  ■  «  .  In  this  case 


x^P{|x^(  >x}  ■  o(  j  y^dP{|x^|  <  y}) 
ly|<x 


and 


lijn  lira  X  f  (s,x)^P{x.  e  b(X)dx}  *  lim  lim  X  f  (s,x)^P{x,  c  b(X)dx} 

e»0  X-^  1  1^  e-*-0  X-M®  I  I  ^ 

lxl<e  lx|<e 

2 

Thus,  noting  u,(2)  *  j  x  dP{jxj^|  <  x}  ,  we  have 

0 

OB 

xj  j  (s,x)^P{exp  {Au}Xj^ €  b(X)dx}du  £ 

^  |xl<e 
/■* 

cX  j  j  (x,x)P{exp  {Au}Xj^e  b(X)dx)du  ^ 


t  I 


X  <e 


-2  /■" 

2cxb  (X)  j  J  {au}}ydydu  * 

-2  j»,Cieb(X)exp{au}  , 

2cXb  (X)  j  j  P{|Xj^|^v}cj  exp  {-2au}dvdu  » 

oo 

c^Xb  ^(X)  J  u2(Cj^eb(X)  exp  {au})  exp  {-2au}du  - 
t 

/oo 

u,(z)z'^e^b^(x)dz  < 

cblAjexp{at} 

c^Xe^U2(eb(X)  exp  {at})(eb(X)  exp  {at))'^  -*■  exp  {-2at}  (3.23j 

(We  have  used  theorems  1  and  la  from  (2.  p.  312314]  which  give  us  the 
following  properties  of  • 

a)  u2(2)  is  a  slowly  varying  function  and 

b)  Xb  ^(X)u-,(b(X))  -►  as  X  -►  ®  } 


Q(s) 


Now  condition  (3.18)  follows  fron  (3.20)  and  (3.23). 

Remark  3.2.  It  follows  from  the  proof  of  theorems  3.2  and  3.3  that 
the  norming  functions  b(X)  and  bj^(n)  for  vrtiich  the  distributions  of 
(x(t,A)-a(A))/b(X)  and  (Xj^+. .  .+Xj^-ajj)/b2^(n)  weakly  converge  can  be 
chosen  eqvial:  b(n)  =  bj^(n)  . 

4.  Limiting  Behavior  of  the  Stationary  Distribution  as  X  -»•  0 
In  this  part,  we  will  suppose  that  the  matrix  A  is  similar  to  the 
diagonal  matrix  A  =  H  ||  ,  vdiere  Xj  ,  1  ^  i  1 »  are  the  eigen¬ 

values  of  A,  i.e. .  A  =  TAT'^  ,  and  T  is  a  nonsingular  matrix  with  real- 
valued  elements.  In  this  case,  Xj^<0,  i^l^n.  We  will  show  that 

the  module  of  x(t,X)  tends  to  zero  with  an  exponential  speed. 

P 

Lemna  4.1.  If  X  -►  0  ,  then  x(X)  -*•  0  .  The  proof  follows  from 
formula  (2.4). 

Let  us  denote  n  *  (rij^ . n^^)  •  t’^Xj^  , 

C  “  (Cl,...,Cjj)  ■=  T'^x(t,X)  , 

Pi  -  Plbi  -  0}  ,  sign  X  *  (sign  Xi,...sign  x^)  if  x  -  (Xi,...Xj^)  , 

Vj  .  . 

For  simplicity  we  consider  only  the  particular  case  Pi^Ojljtij^n. 

Theorem  4.1.  If  Pi  ■  0  for  all  i,  1  i  ^  n  ,  then  the  distribution 
of 

(sign  ;  ,  ICjl 

weakly  converges  as  X  -*•  0  to  the  distribution  of 

(sign  n  ,a, . . .a) 


where  a  has  the  uniform  distribution  on  the  interval  (0,1)  and  does  not 
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depend  on  n. 

Proof.  The  distribution  of  x(t,X)  coincides  with  the  distribution 
of  the  vector 

C  =  J  exp  {Au}dz(u)  (4.1J 

0 

We  can  suppose  that  the  process  z(t)  is  determined  by  the  values  of 
its  juin)s  Xj^,  X2,...  and  by* the  lengths  of  the  intervals  between  junps 
X  X  where  all  x^  and  are  independent, 

P{t^>x}  =  exp  {-x},  X  >  0  . 

Thus,  the  formula  (4.1)  inplies 

C  =  exp  {X  ^At2}xj^  +  exp  (X'^ACtj^  +  T2) }x2  +...  =  exp  {X*^ATj^}(Xj^  +  C^)  , 

where  ,  and  ^  are  independent  and  the  distributions  of  5  and 

coincide.  It  follows  from  (4.2)  that 

T'^5  =  exp  {x'^Atj^}T'^(x^  +  C^) 

or 

K  »  exp  (X  ^Atj^Xn  +  K^)  where  <  =  . 

ic^  -  and  the  distributions  of  ic  and  c  coincide  (4.3) 

According  to  lemma  4.1  0  and  we  have  from  (4.3) 

1  P 

sign  <  ■  sign  (n+x  )  sign  n  , 

kjl  "  |exp  (X'^X^T^Xn^  ^kJ)!  ^  exp 

Since  exp  has  uniform  distribution  on  (0,1)  and  does  not 

depend  on  n,  the  statement  of  the  theorem  easily  follows  from  the 


(4.2) 
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well -known  properties  of  convergence  in  probability  and  weak  convergence. 

5.  Applications  to  a  Queueing  System  with  Changeable  Service  Rate 

In  this  part,  we  consider  the  case  n  =  1  ,  Xq  ^  0  ,  ^  0  .  We 

write  equation  (1.1)  in  the  form 

dx(t)  =  -ux(t)dt  +  dz(t)  (5,1) 

where  y  >  0  . 

Equation  (5.1)  is  connected  with  a  following  queueing  system.  Input 
flow  is  a  Poisson  flow  with  parameter  X.  To  serve  the  n’th  customer,  the 
server  has  to  spend  units  of  work.  If  x(t)  is  the  amount  of  work 
necessary  to  serve  the  customers  present  in  the  system  at  the  moment  t,  then 
service  rate  equals  to  yx(t) .  We  will  investigate  the  distribution  of  the 
virtual  waiting  tune  0  subject  to  the  condition  that  system  has  stationary 
distribution.  (It  is  assumed  that  system  has  FIFO  service  discipline.) 

Let  y(t)  be  the  amount  of  work  performed  by  the  moment  t.  Then 

y(0)  =  0  , 

dy(t)  =  ux(t)dt  (5.2) 


Thus 


dy(t)  *  dz(t)  -  dx(t) 
and 

nt)  .  i(t)  -  x(t)  .  .  Xod-e-"') 

0 

where  Xq  has  the  distribution  defined  by  (2.4). 
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If  we  denote  by  x  the  amount  of  work  necessary  to  serve  the  considered 
customer  then 

P{e>t}  =  P{y(t)  <  Xg/x(0+)=  Xq  +  x}  = 

P{  (l-e‘''^^*^bdz(u)  +  (XQi-x)  <  x^}  = 

0 

P{  (l-e‘^^^'^^)dz(u)  +  x(l-e'^^)  <  e'^’^Xp}  (5.3) 

0 

Theorem  5.1.  If  Ex  =  iaj^<®,  then  y0  1  as  X  -*•  »  . 

Proof.  The  proof  follows  from  the  relations 

X  (l-e'^f^^'^^JdzCu)  (l-e'^‘^^'“bm^du  = 

0  0 

mj  -  m^y‘^Cl-e‘^^)  , 

e^^<iziu)  e^\du  =  y'V^^  , 

Jbo 

X*^x(l-e'^^)  ^0  as  X  -*•  “  . 

Theorem  5.2.  If  the  distribution  of  x  belongs  to  the  domain  of  attrac¬ 
tion  of  a  stable  law  with  exponent  a,  1  <  a  <  2  ,  then  there  exists  a 
regularly  function  f(X)  with  exponent  1-a’^  such  that  the  distribution 
of  (0-l/v)f(X)  weakly  converges  as  X  -*■  ®  to  the  stable  law  with  exponent 
a  and  c.f. 

h(s)  «  exp  {(e^’^“^^(ow)'^{mj^e)'“  + 

-1 

/  e-"“/^l-e-““)“<iu)s'‘)  ,  s  >  0  . 

0 


Proof  is  analogous  and  follows  from  the  theorems  3.2,  3.3  and  repre¬ 
sentation  (5.3). 

Theorem  5.3.  If  the  distribution  of  x  belongs  to  the  domain  of 
attraction  of  a  stable  law  with  exponent  a,  0  <  a  <  1  ,  then 

lim  P{0  >t}  =  PCCZj^/Zp"  >  af  (e^''^-e^^“)du}  , 

0 

vdiere  Z^  and  Z2  are  positive  independent  identically  distributed  random 
variables  and  have  stable  distribution  with  exponent  a. 

Proof  is  analogous  and  follows  from  the  theorem  3.2  and  representation 
(5.3). 

Theorem  5.4.  lim  P(0^^^  <  t}  =  t  ,  if  t  S  (0,1)  . 

Proof  follows  from  the  representation  (5.3)  and  theorem  4.1. 
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